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Abstract A certain vector-tensor theory is revisited. 
Our attention is focused on cosmology. Against pre- 
vious suggestions based on preliminary studies, it is 
shown that, if the energy density of the vector field 
is large enough to play the role of the dark energy 
and its fluctuations are negligible, the theory is not si- 
multaneously compatible with current observations on: 
supernovae, the cosmic microwave background (CMB) 
anisotropy, and the power spectrum of the energy den- 
sity fluctuations. However, for small enough energy 
densities of the vector field and no scalar fluctuations, 
the theory becomes compatible with all the above ob- 
servations and, moreover, it leads to an interesting evo- 
lution of the so-called vector cosmological modes. This 
evolution appears to be different from that of general 
relativity, and the difference might be useful to explain 
the anomalies in the low order CMB multipoles. 

Keywords cosmology:theory~large-scale structure of 
universe-modified theories of gravity 



1 Introduction 



2008ll . 



In previous papers ([Morales fc Saea 12007 . 
number of effects produced by superhorizon cosmolog- 
ical vector modes were discussed in the framework of 
General Relativity (GR); in particular, it was proved 
that these modes may explain the anomalies in the first 
multipoles of the cosmic microwave background (CMB) 
temperature distribution. In GR, vector modes decay 
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during the matter dominated era and, consequently, it 
is not easy to explain their presence at redshifts close 
to z ~ 1100 (matter-radiation decoupling); however, 
in vector-tensor (VT) theories, the evolution of vector 
modes is expected to be different and, in some cases, 
this evolution might be appropriate to explain CMB 
anomalies. It occurs in the case of the VT theory 
studied in this paper; in fact , as it is proved below 
(see also Dale. Morales fc Saez ( 20091) ), for appropriate 
small values of the vector field energy density without 
scalar perturbations, the theory is compatible with cur- 
rent observations and, moreover, superhorizon vector 
modes grow during the matter dominated era. 

It is well known that the so-called concordance model 
simultaneously explains most of the current cosmolog- 
ical observations. By this reason, this cosmological 
model is considered -along the paper- as the preferred 
model in the framework of GR. In the concordance 
model, the universe is quasi-flat and the initial fluc- 
tuations have an inflationary origin; in this situation, 
scalar perturbations are fully dominant, whereas the 
effects due to gravitational waves (vector modes) are 
expected to be small (negligible). Moreover, scalar 
perturbations are purely adiabatic, their distribution 
is Gaussian, and the initial power spectrum is simi- 
lar to a Harrison-Zel'dovich one. There are a set of 
observations to be simultaneously explained by a con- 
cordance model; e .g.: (i) WMAP observations of the 
CMB anisotropics (Komatsu. et aDbOfd Jarosik. et al 



2OIOI: iLarson. et allbolOl) . (h) other CMB anisotropy 
observations -involving greater ^-values- performed 
from ground and ballo on experiments; amon g them 



( Fowler et al. 


2OIOI) and SPT (Lueker et al. 


2OO9I). 


(iii) high-redshift la supernovae (SNe la) observations 


(Astier et al. 2006: Riess et al. 


200fil: 


Kowalski et al. 



20081 : iFolev et al.ll2009l) . Ti 



powe r spectra measure - 
ments based on galaxy surveys (|Reid et all boOOl) . 
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(v) detection of baryonic acoustic oscillations (BAO) 



(Eisenstein et al.lboOSl Takahasliil2009l 


Martinez, et al. 


20091 Percival et alJ 2009). (vi) measi 


irements of the 
id (vii) primor- 


Hubble constant (Suvu et al. 


2010|), ai 



dial deuterium abundance observations (jPettini et al, 
I2OO8I) . In the LAMBDA archive (WMAP seven years 
data), various models -including different effect as 
lensing, Sunyaev-Zel'dovich, and so on- are fitted 
to different sets of observational data including al- 
ways the proper WMAP7 data; e.g., for a ACDM 
model including lensing and Sunyaev-Zel'dovich ef- 
fects, the best fitting to WMAP7 plus BAO and Hub- 
ble constant nieasu rements corresponds to (see also 
Jarosik. et all |2010l )): (1) a reduced Hubble constant 



h = IQ-^Ho = 0.704t°:°^i^ (where Ha is the Hubble 
constant in units of km s"-^ Mpc^^); (2) density pa- 
rameters nb = 0.0456 ± 0.0016, r2c = 0.227 ± 0.014, 
and fl\ ~ 0.728lg'Qi[g for the baryonic, dark matter, 
and vacuum energy, respectively (the matter density 
parameter is then f2m = fJfc -I- ~ 0.2726); (3) a total 
density parameter fi^ = 1.0023lg;Qg54, (4) a parameter 
(78 = 0.809 ± 0.024 normalizing the power spectrum of 
the scalar energy perturbations, (5) a running scalar 
spectral index = 0.963 ± 0.012 with dns/d{ln k) = 
-0.022 ± 0.020 for fc = fcg = 0.002 Mpc'^, (6) an 
equation of state for dark energy of the form p = Wp 
with W = -0.980 ± 0.053, and (7) an optical depth 
T — 0.087 ± 0.014 which characterizes the reionization. 
There are fifty eight different fittings in the LAMBDA 
archive. All of them give similar CMB angular power 
spectra to fit WMAP7 data. By this reason, results 
in next sections appear to be almost independent of 
the particular choice of the fitting parameters, provided 
that they are close enough to the central values of some 
LAMBDA fitting. Our choice of these parameters is 
done in section [3] and, then, the resulting concordance 
model is assumed along the paper as a model of refer- 
ence leading to a good fitting with current observations. 
Predictions from models based on the VT theory un- 
der consideration will be compared with observational 
data. Obviously, only predictions similar enough to 
those of the reference model may be compatible with 
current observations. 

This paper is structured as follows: The VT the- 
ory is described in Section [2l predictions of the concor- 
dance model are compared with observations in section 
[3] These comparisons are extended to VT cosmology in 
section m Vector perturbations are studied in section [S] 
and, finally, section [5] is a general discussion and a sum- 
mary of the main conclusions. Let us finish this section 
fixing some notation criteria. Latin (Greek) indexes 
run from 1 to 3 (0 to 3). The gravitational constant 
and the scale factor are denoted G and a, respectively. 



Units are chosen in such a way that the speed of light 
is c = 1. 



2 The theory and its basic cosmological 
equations 

In vector-tensor theories, there are two fields, the metric 
g^v and a four- vector A>^. Various of these theories 



have been proposed (see IWilll (|l993l ). IWilll (|2006[ ) and 
references cited there). We a re intereste d in one of the 
theories based on the action ( Will 19931) : 



/ = il6TTG)-^ J{R + ujA^A>'R + r]Rf,^A''A'' - 

eF^.F^'' +-fV,A^V''A^ + L^)^d^x (1) 

where w, rj, e, and 7 are arbitrary parameters, R, R^u, 
g, and Lm are the scalar curvature, the Ricci tensor, 
the determinant of the g^^i, matrix, and the matter La- 
grangian, respectively. The symbol V stands for the 
covariant derivative, and F^i, = V^A^ — S/^A^. In 
action ([TJ, it is implicitly assumed that the coupling 
between the matter fields and A,, is n egligible. A dif- 
ferent action is given in (j2006[ ). It involves a 
new term of the form X{A^A^^ + 1), where A is a La- 
grange multiplier. With the help of this term, the vec- 
tor A^ is constrained to be timelike with unit norm. 
From this last action, the field equations of the so- 
called Einstein- Aether constrained theories can be eas- 
ily obtained. The applicatio ns of these theories to Cos- 
molog y are discussed, e.g., in lZuntz. Ferreira fc Zlosnik 
( 2008h. A mass t e rm of the form m\A^A^ is used by 
Bohmer fc Harkol (|2007t ) to expla in cosmic acceleration 
with a massive vector field. The same acceleration was 
explaine d in the framework of the the o ry stu died in this 
paper bv iBeltran Jimenez fc Marotd ([20081). Recently, 
other theories involving vector fields have been also 
apphe d to cosmology (see, e.g 
( 20071 ): iMoffatI (|2006h 'l. 



Tartaglia fc Radicella 



Whatever the action parameters may be, only the 
PPN parameters 7, /3, ai and a2 may be different from 
those of GR in the theories based on action ([1]). More- 
over, in these theories there is a ya.rying effective grav- 

On account 



Will (199 



itational constant Ge// (see 
of current solar system observations, we are only in- 
terested in theories having 7 = /? = 1 and Geff — G. 
This condition is hereafter our PPN viability condition. 
For w = £ = 0, there are only two theories compati- 
ble with the PPN viability condition, the first (second) 
of these theories corresponds to = —7 (77 = 7). In 
the second case, we have verified that vector cosmo- 
logical modes evolve as in GR; hence, on account of 
our motivations (see section [T]) , we are not interested 
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in this theory but in a theor y with 77 = —7. From 
the PPN parameters given in IWilll (1993), it is easily 
seen that any theory with the action parameters cj ^ 0, 
arbitrary 77, 7 = 77, and e = — w + ri/2 is compatible 
with the PPN viability condition. Finally, the same oc- 
curs for w = 0, arbitrary e and 7 = 77 and also for 
w = 0, arbitrary e, and 77 — 4e — 7. Although all 
these theories deserve attention, their study is out of 
the scope of this paper. Here, we only analyze the the- 
ory with 77 = —7 = Stt G* and a; = £ = 0, whose fiel d 



equations were given in iDale. Morales fc SaezI (|2009l ). 
Various reasons have motivated the choice of this the- 
ory: (1) it is the simplest theory compatible with the 
PPN viability condition, (2) it does not involve new 
dimensional constant s, (3) it has been already stud- 
ied in various papers (iBeltran Jimenez fc Marotollioosi 
Beltran Jimenez. Lazkoz fc M aroto 2009) , where it was 
proved that -in this theory- there are no classical and 
quantum instabilities, and also that the so-called cosmic 
coincidence problem is avoided, and (4) deeper study is 
necessary to support its reliability. 

Scalar and vector A'^-perturbations of cosmological 
interest should arise inside the effective horizon dur- 
ing inflation; nevertheless, the VT theory does not de- 
fine the inflationary period, which is likely produced by 
some scalar field. The coupling between this scalar in- 
flaton and the vector field is basic to estimate the 
initial inflationary ^^-perturbations. This coupling is 
to be chosen among many possibilities. After reheating, 
the scalar perturbations evolve coupled to the remain- 
ing scalar modes (associated to the metric, the energy 
densities, and so on) and their evolution is very compli- 
cated. A code similar to CAME, but including the new 
scalar ^'^-modes, would be necessary. This study and 
the full analysis of possible inflationary models is be- 
yond the scope of this paper. On account of these facts, 
we first complete the study of the selected VT theory 
in the case of negligible ^^-perturbations (see previ- 
ous papers by Beltran Jimenez fc Marotol (|2008) and 
Beltran Jimenez. Lazkoz fc Marotd (j2009f )) and, after- 
wards, we study the evolution of possible vector pertur- 
bations after reheating. Encouraging results from this 
study (see sections [5] and [6]) suggest future work (see 
section [6]) in the framework of both the chosen VT the- 
ory and other VT theories compatible with the PPN 
viability condition. 

From the field equations (|Dale. Morales fc Saez 

one easily finds the equations describing a ho- 



mogeneous and isotropic universe, in which the line 

element is 



-dr' 



1 -Xr2 



2 • 2 
r sm 



and the vector field has the covariant components 
(Ao(t), 0, 0, 0). Constant K takes on the values +1, —1 
and in closed, open and flat universes, respectively, 
and r is the conformal time. The resulting equations 
are: 



a 



3— = SttGo^ -I- + 



-2- + — = S^Ga^ K +Pa . 
where, 



Va 



a 



4d 



2a4^° 



a" 



2 2 



^0 ' 



Pk =-^Pk = -SK/SirGa^ , 



(3) 

(4) 
(5) 

(6) 

(7) 
(8) 



(2) 



the dot stands for a derivative with respect to the time 
r, and and are the energy density and pressure of 
the cosmological fluid, respectively. This fluid contains 
baryons, radiation, and cold dark matter (CDM). We 
see that the field equations of the theory couple the 
variables Aq and a describing the background. The 
density parameters due to baryons, radiation, CDM, 
the vector field and the curvature term are fib, fir, ^c, 
ri^ , and , respectively. They satisfy the relation 
fib + + flc + + = 1- In the flat case one has 
ftjf — and ~ rtb + ftr + fic + ^A whereas in 
the closed (open) case the curvature density parameter 
is negative (positive) and the inequality — flb-\-^r + 
flc + flji > 1 {ft < I) is satisfied. 

The free parameters of the theory are assumed to 
be the Hubble constant and the density parameters fi^ 
and flj^ . In the fiat case, the present value of the scale 
factor, flo, is arbitrary (we assume oq = 1); whereas in 
nonflat models this value is gq = 7Jq"^|J7^^ j^^/^^ por a 
given choice of the free parameters, Eqs. ([3]) - ([5]) may 
be numerically solved to get functions a — a(T) and 
Aq = Ao{t). Finally, by using Eqs. dl]) - 0, the equa- 
tion of state for the dark A^-energy; namely, the ratio 
W{t) — Pa{t)I PAi"^) be easily calculated. We can 
then state that, in the absence of A^-perturbations, 
the VT theory is cosmologically equivalent to GR plus 
dark energy with the computed W{t) ratio. This vari- 
able W{t) may be e asily used for numerical calculation s 
with CAMB code |Lewis. Challinor fc Lasenbvl I2OO0I) . 
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This code allows us to calculate all the cosmological 
spectra, including those associated to the CMB polar- 
ization. All these spectra are influenced by the vector 
field of the VT theory. The effects of the vector field 
would be identical to those produced, in GR, by dark 
energy with the associated equation of state {W{t) ra- 
tio computed in the VT theory). 



3 Supernovae, CMB and Power Spectrum 
Observations 

The concordance model is now defined and analyzed in 
detail. The contents of this section are used along the 
paper to discuss the viability of the vector-tensor theory 
under consideration. Any admissible theory must ex- 
plain the same observations as the concordance model 
with comparable accuracy. Hereafter, this model is de- 
fined by the following set of parameters, which are close 
to the central values for the fitting considered in section 
II: h=0.704, = 0.0461, ^ 0.2265, = -0.0012, 
r = 0.087, W = -1, Us = 0.96, dns/d{ln k) = -0.016 
and as = 0.8065. These parameters are used to get the 
angular correlations of the CMB and the power spec- 
trum P{k) of the matter fluctuations. Calculations are 
performed by using CAME (see section [5]) . The rela- 
tion between the distance modulus fi — m — M and 
the redshift z of the la supernovae only depends on the 
parameters h, fi™ and fifc. In other words, this rela- 
tion only depends on some parameters describing the 
background universe, whereas the parameters associ- 
ated to perturbations are irrelevant. Let us now study 
the agreement of this model with data from various ob- 
servations. 

For the concordance model (h=0.704, r2„i — 0.2726 
and fifc = —0.0012), the relation /x — fj.{z) is repre- 
sented in the solid line of Fig. [TJ In the same Figure we 
exh ibit 156 observatio nal data, from lAstier et al. I (l2006l) 
and lRiess et al.l (|2006l ) , which cover a wide range of red- 
shifts from z = 0.015 to z = 1.755. Error bars do not 
include the uncertainty in the absolute magnitude of 
the SNe la, which is typically 0.1 5 mag, we t h en ca l- 
culate the following function (see Astier et al. ( 20061) 1 



E 



(9) 



where /^obs ilJ-the) is the observational (theoretical) /i 
value, aint accounts for the intrinsic dispersion of SNe 
la and a{fi) is the error represented in Fig. [TJ For 
Cint = 0.15 mag, the ^^^'^^ is 132.45 and the p- 
value is 0.91. This is a large value suggesting that 



the chosen parameters explain very well the SNe la ob- 
servations. Of course, these parameters are very sim- 
ilar to those of section [1] which are compatible -by 
construction- with WMAP7, BAO, and Hq observa- 
tions. 

Let us now study the CMB angular power spec- 
trum for the concordance model. CAMB calculations 
are performed including lensing. The resulting spec- 
trum is given in the solid line of Fig. [2l where we 
also present the binned WMAP seven years measure- 
ments with the corresponding error bars, which can 
be found in LAMBDA (Legacy Archive for Microwave 
Background Data Analysis). The comparison between 
observations and theoretical predictions is performed 
by using the function xl^^ = Y^Mnsi^obs - Athe)^/a^, 
where A = £{£+l)Ce/2Tr and Aobs i^the) stands for the 
observational (theoretical) A-values, respectively. The 
dispersion of A inside each £-h'm is measured by a^. 
From the 45 bins considered by the WMAP team, the 
resulting Xcmb 36.81 and the corresponding p-value 
is 0.80. This high p-value strongly suggest that the con- 
cordance model fits very well WMAP7 observations. 

The spectra corresponding to the electric part of the 
CMB polarization and the cross correlation TE have 
been also obtained, but they are not necessary for the 
discussion presented in this paper. 

The matter power spectrum estimated with CAMB 
is shown in the solid line of Fig. [51 to gether with ob- 
servational data ([Eisenstein et al.ll2005[ ) from the sloan 
digital sky survey (SDSS). Baryonic acoustic oscilla- 
tions (BAOs) would produce the deviation between 
the solid (CAMB spectrum) and the dashed lines of 
the top panel of Fig. [4] The dashed line shows the 
spectrum in the absenc e of baryon acoustic osci llations 
|Hu fc EisensteinI Il998t lEisenstein fc Hul Il998l) . Par- 
ticularly relevant is the deviation appearing in the k/h 
interval (0.05, 0.1). Recent obse rvations seems to sup- 
port the existence of this feature ( Eisenstein et ahlEoOSt 
Martinez, et al.l 120091 : IPercival et al.ll2009t ). which has 



been taken into account (by the WMAP team) to find 
the parameters of section [T] 

As it is shown in Fig. [Sj the observed and theoretical 
power spectra fit very well for large spatial scales with 
k < 0.05/i Mpc~^, which are clearly linear. For smaller 
scales, the error bars of the observational data are above 
the solid line (predictions based on the concordance 
model). The theoretical spectrum is due to perturba- 
tions in both dark and baryonic matter, whereas the 
power observed in galaxy surveys corresponds to fluc- 
tuations in the baryonic component alone. These two 
spectra would be only identical in the absence of any 
bias between dark and baryonic energy densities at the 
corresponding spatial scales. Hence, if there are no un- 
known systematic errors in the observations. Fig. [3] 
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Table 1 Parameter configurations in GR and VT cosmologies. 



Model 




h 










T 


0-8 


Us 


dns/d{ln k) 






GRl 


.2726 


.704 


.0228 


.1123 


.7286 




.087 


.8065 


.96 


-.016 


-.0012 


36.81 


VTl 


.25 


.82 


.0223 


.1458 


0. 


.7512 


.087 


.8065 


.96 


-.016 


-.0012 


15917. 


VT2 


.25 


.82 


.0223 


.1458 


0. 


.7512 


.000 


.95 


.96 


-.016 


-.0012 


2467. 


VT3 


.25 


.82 


.0223 


.1458 


0. 


.7512 


.000 


.95 


.92 


-.016 


-.0012 


1473. 


VT4 


.25 


.82 


.0223 


.1458 


0. 


.7512 


.000 


.95 


.84 


-.016 


-.0012 


432.8 


VT5 


.25 


.82 


.0223 


.1458 


0. 


.7512 


.000 


.95 


.84 


-.05 


-.0012 


288.4 


VT6 


.25 


.82 


.0223 


.1458 


0. 


.7512 


.000 


.95 


.84 


-.09 


-.0012 


180.9 


VT7 


.25 


.82 


.0223 


.1458 


0. 


.7700 


.000 


.95 


.84 


-.09 


-.02 


1057.9 


VT8 


.25 


.82 


.0223 


.1458 


0. 


.7300 


.000 


.95 


.84 


-.09 


.02 


355.0 


VT9 


.25 


.82 


.0223 


.1458 


0. 


.7400 


.000 


.942 


.84 


-.09 


.01 


83.28 


VTIO 


.36 


.73 


.0223 


.1695 


0. 


.6412 


.087 


.8065 


.96 


-.016 


-.0012 


17725. 


VTll 


.36 


.73 


.0223 


.1695 


0. 


.6412 


.000 


.95 


.96 


-.016 


-.0012 


3802. 


VT12 


.36 


.73 


.0223 


.1695 


0. 


.6150 


.000 


.963 


.84 


-.13 


.025 


86.62 


VT13 


.48 


.67 


.0223 


.1932 


0. 


.5212 


.087 


.8065 


.96 


-.016 


-.0012 


19683. 


VT14 


.48 


.67 


.0223 


.1932 


0. 


.5212 


.000 


.95 


.96 


-.016 


-.0012 


5412. 


VT15 


.48 


.67 


.0223 


.1932 


0. 


.4800 


.000 


.993 


.84 


-.16 


.04 


136.9 


VT16 


.2726 


.704 


.0228 


.1123 


.7276 


.0010 


.087 


.8065 


.96 


-.016 


-.0012 





suggests a small bias, for k > 0.05/i Mpc~^, whose 
origin has not been explained: We might speculate 
with some kind of nonlinear effect, with a primordial 
bias depending on the spatial scales (which vanishes 
for k < 0.05h Mpc~^), and so on. 

Another characteristic of the matter power spectrum 
is the CTg value, which is essentially related to the power 
at scales smaller than 8h~^ Mpc; namely, to the power 
corresponding to wavenumbers k > 0.12bh Mpc. Var- 
ious estimates of erg have been reported in the techni- 
cal literature. After seven years of CMB observations, 
WMAP team has obtained the value erg ~ 0.8. More- 
over, recent measurements of the CMB a ngular power 



spect r um, at very small angular scales (jFowler et al 



2OI0I : iLueker et all l2009t ). suggest CTg values smaller 
than ^ 0.86. Similar bounds have been obtained 
from the ana lysis of ROSAT X-ray cluster data (see 
Fowler et al.l ([2OIO) and references cited therein). Fi- 
nally, the CTg value corresponding to the SDSS data of 
Fig. [3] is close to 0.84. On account of all these re- 
sults, we hereafter accept the observational constraint 
(Tg < 0.9 for dark plus baryonic matter. 

For a given theoretical model, the resulting p- values 
depend on the chosen observational data, but the data 
used along this section have been selected among the 
most accurate current observations (WMAP for the 
CMB, SDSS for the galaxy distribution, and the su- 
pernova legacy survey and Hubble observations for su- 
pernovae) and, consequently, we can be confident with 
our results and use the same data to study the compat- 
ibility between theoretical predictions and observations 
in the framework of the VT theory under consideration. 
It is done in next section. 



4 Results 

We begin with the study of the apparent supernovae 
dimming in the framework of the VT theory described 
in section [2j In this theory and also in any theory with 
cosmological models based on the Robertson- Walker 
background metric, CMB observations ensure that the 
universe is almost flat. It is due to the fact that the 
first peak of the CMB angular power spectrum appears 
located at £ ~ 210. On account of these facts, the value 
of the curvature density parameter has been chosen to 
be fife = —0.0012 as in our version of the concordance 
model. In the VT theory, CAMB numerical calcula- 
tions can be easily performed in the presence of this 
small curvature and, consequently, it has been main- 
tained all along this section. Nevertheless, it has been 
verified that, with the chosen ri^-value, results are al- 
most indistinguishable from those corresponding to a 
strictly flat universe (see first paragraph of section [5]) . 
Parameters h and flm have been varied. Only parame- 
ters f2fe, h, and f2„i are relevant in supernova studies. 

For each pair {h, ftm), function Xcmb ^^'^ 
sociated p-value have been calculated. The same ob- 
servational data as in section [3] have been used in the 
computations. From a purely statistical point of view, 
models leading to p < 0.05 are usually ruled out. In 
Fig. [5] you can see a 3D representation involving quan- 
tities h, Qm, and lOOp (p-value %). In this Figure, we 
see that the p-value is larger than 0.05 in a bounded 
region of the (h, flm) plane. This region is exhibited in 
Fig. [6l where the lOOp-values in the {h, plane are 
represented by using a color scale. For any small fifc pa- 
rameter allowed by CMB observations, the shape and 
size of this region are very similar. According to Fig. 
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[71 condition p < 0.05 rules out f2„i values smaller than 
^ 0.24. We also see that fim-values a little greater than 
0.48 are not discarded by the same condition; however, 
they will be discarded from our analysis of the CMB 
angular power spectrum (see section!?]). 

Let us now study the CMB angular power spectrum 
for all the pairs {h, flm) explaining supernovae observa- 
tions. We have studied forthy two {flm, h) pairs corre- 
sponding to seven h- values and six f2m-values uniformly 
spaced in the intervals [0.6, 0.9], and [0.25, 0.55], respec- 
tively. These pairs cover the region of the h-flm plane 
allowed by SNIa observations for any small ri^-value 
compatible with the first peak location in the CMB an- 
gular power spectrum. We cover a region a little greater 
than that described in Figs. [3 El and [71 Each (fi™, h) 
pair is studied by using the same method, which is now 
described in detail for three selected additional pairs 
located in the region of Fig. [6l In the top left panel of 
Fig. [H we show Ci quantities for the pair {h — 0.82, 
= 0.25), which poorly explains supernovae observa- 
tions with p ~ 0.07. The middle left panel corresponds 
to h = 0.73 and r2„i — 0.36 and supernova observations 
are well explained with p ~ 0.82 and, finally, in the bot- 
tom left panel we have represented the case h = 0.67 
and rim = 0.48 having p ~ 0.31. These three pairs cor- 
respond to the maxima of the curves presented in Fig. 
[7] for = 0.25, n,n = 0.36, and = 0.48. 

Let us begin with the pair h = 0.82, flm — 0.25, 
which is not a good choice from th e point of view of s u- 
pernovae and h measurements (see Suvu et al. ( 2010l )): 
however, it appears to be the best choice to deal with 
the CMB (see below). By this reason, this pair plays a 
central role in our discussion on CMB anisotropy. The 
dotted line of the top left panel of Fig. [51 is the an- 
gular power spectrum for the parameter configuration 
VTl of Table [H In this configuration, the value of 
flfjh^ is assumed to be 0.0223, which is the maximum 



possible values are considered in the last paragraph of 
this section. 

We can now modify some of the remaining parame- 
ters to look for a good set of Ci coefficients. Actually, 
the following parameters may be changed; r, erg, fi^, 
Ug, and dus/ d{ln k). Changes in the parameters r and 
(Tg may produce a large magnification of the angular 
power spectrum amplitude. Since this amplitude in- 
creases as CTg (r) increases (decreases), the values r = 
and (Tg = 0.95 lead to a relevant magnification, which 
is not realistic as a result of the following facts: (i) 
the inequality erg < 0.9 should be satisfied (see sec- 
tion IS]) and, (ii) reionization at a redshift z ~ 6 is 
proved by measure ments of the Gunn-peterson effect 
( Becker et al ] l200lh . which implies r > 0. For the cho- 
sen new extreme values of r and erg, plus the values 
Us = 0.96 and dns/d(ln k) — —0.0016 of the concor- 
dance model (parameter configuration VT2 of Table [Ij, 
the Ci quantities have been calculated, they are shown 
in the dashed line of the top left panel of Fig. [8] We 
see that, in spite of our forced choice of r and erg, the 
resulting angular power spectrum remains smaller than 
that of the concordance model (solid line). In this case, 
we have found Xcmb ~ 2467. and p ~ 0. 

The same study has been performed for the pairs 
{h = 0.73, Urn = 0.36) and {h = 0.67, = 0.48). 
Results are presented in the left middle and left bot- 
tom panels of Fig. [51 In the middle left (bottom left) 
panel the dotted line is found for the configurations 
VTIO (VT13) of Table [H whereas, the dashed hues 
are obtained for the configuration VTll (VT14). As 
in the top left panel, the dashed lines correspond to 
r = and cg = 0.95. For these lines, we have found 
ixL.s = 3802., p ^ 0.) and (x^,,^ = 5412., p ^ 0.) 
for the middle left {h — 0.73, fim — 0.36) and the bot- 
tom left {h — 0.67, n„i = 0.48) panels, respectively. 
These -values are larger than that corresponding 



possible valu e at 6 8 % confidence level according to to h — 0.82 and ilm — 0.25 (x? 



2467.); hence, the 



Pettini et al.l (|2008f ). the ftc value is then fixed by the 
relation flm = ilb + — 0.25 and, the remaining pa- 
rameters are taken identical to those of the concordance 
model. In the same panel, and also in any left panel of 
Figs. [5HTU1 the solid line shows the angular power spec- 
trum of the concordance model (GRl entry of Table[T]). 
The same is valid for any right panel, where the solid 
line is the matter power spectrum of the case GRl. 
We easily see that the chosen parameters lead to very 
small Ct quantities (dotted line) which cannot explain 
WMAP observations. Accordingly, for these parame- 
ters and the observational data of Fig. [21 we have ob- 
tained Xcmb ~ 15917. and p ~ 0. . In configurations 
VT1-VT15 of Table [H th e parameter f^ bfe^ takes on its 
maximum possible value iPettini et al. f200&l . Other 



best situation is found for this last pair. The same oc- 
curs for any pair compatible with SNe la observations. 

After these considerations we continue with the 
study of the best pair. The question is: Could we fit 
the observations varying $7^, Us, and dns/d{ln k) in the 
parameter configuration VT2? We begin with parame- 
ter Us, whose value, in VT2, is = 0.96. In the new 
configurations VT3 and VT4 we have taken = 0.92 
and ria — 0.84, respectively. The remaining parameters 
have not been altered (see Table [1]). The resulting an- 
gular power spectra are exhibited in the top left panel 
of Fig. [9l The dashed (dotted) line corresponds to 
n, = 0.92 {us = 0.84). For the VT3 (VT4) model we 
have found xl,,^ = 1473. (x^„^ = 432.8) and p ~ 0. 
[p ~ 0.). Accordingly, we see that the dashed line is 
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well below the solid line. The dotted line approaches 
rather well the first peak, but it does not fit the low I 
multipoles and the remaining peaks of the solid line. A 
good fitting based on parameter Us is not possible for 
any I < 1150 (the largest £ value in the WMAP seven 
years binned data). 

It is also remarkable that the value — 0.84, which 
leads to a good fit of the first peak, is too small to 
be compatible wi th standard infla- tionary models, for 
which. Us > 0.94 ( Page et al. 200?!). Rather ex otic in- 
flationary models ( Mukhanov fc Vikman 20061 ). whose 
Lagrangians involve nonlinear functions of the inflaton 
kinetic energy might be a rare exception. 

Let us now change the value of dns/d{ln k) in config- 
uration VT4 of Table[T] Two new values of dns/d{ln k): 
—0.05 and —0.09 have been assumed in configuration 
VT5 and VT6, respectively. The angular power spec- 
trum of these two new configurations are shown in the 
left middle panel of Fig. 13 The dashed (dotted-dashed) 
line corresponds to the value —0.05 (—0.09). The dot- 
ted line is identical to the corresponding line of the top 
left panel. Comparison with WMAP binned data gives 
xL.B = 288.4 (xL,, = 180.9) for the value -0.05 
(—0.09) and p ~ 0. in both cases. 

As a last step, the parameter ilk has been slightly 
changed in configuration VT6. The new values -0.02 
(closed universe) and 0.02 (open universe) have been 
assumed in the configurations VT7 and VT8 of Table 
m respectively. Results are presented in the bottom left 
panel of Fig. [S] The dotted (dashed) line corresponds 
to the closed (open) case. In the closed (open) case, the 
peaks are lower (higher) than in the concordance model 
and they are shifted to left (right). The low i multipoles 
are almost independent of these changes in ilk ■ Since 
too large peak shifts are inadmissible, parameter flk is 
constrained to be smaller than a few hundredths. 

Finally, we have varied parameters as and flk in 
configuration VT8 to optimize the fitting to the ob- 
servational data. The resulting values as — 0.942 and 
fifc = 0.01 (slightly open case) define the VT9 con- 
figuration. The angular power spectrum of this last 
configuration is displayed in the dotted line of the top 
left panel of Fig. (TUl In cases (h = 0.73, flm = 0.36) 
and (h = 0.67, ilm = 0.48), the best fittings have been 
found for {as = 0.963, flk = 0.025) and {as = 0.993, 
fifc = 0.04), respectively. The resulting spectra are dis- 
played in the dotted line of the middle left {h — 0.73, 
rim = 0.36) and bottom left {h = 0.67, Qm = 0.48) 
panels of Fig. [TUl In the top left panel {h = 0.82, 
= 0.25), we see a rather good fitting to the ob- 
servational data excepting the first ten multipoles. As 
a result of this discrepancy one has Xcmb ~ 83.28 and 
p ~ 4.52 X lO^'*. In the middle left and bottom left pan- 
els the fittings are slightly worse in the regions of both 



the first multipoles and the third peak. Accordingly, we 
have found X^^mb ~ 86.62 for the middle left panel and 
= 136.9 for the bottom left one. All these 

^CMB ^CMB 

values are much greater than 36.81, which is the value 
corresponding to the concordance model. For pairs {h, 
flm) with Qm > 0.48, the fittings are worse. 

The matter power spectra of all the configurations 
considered in the above discussion have been presented, 
with the same type of line, in the corresponding right 
panels of Figs. [H [H and [TOl From these figures, it 
follows that the matter power spectrum and the CMB 
angular power spectrum do not tend to the spectra of 
the concordance model at the same time. The VT con- 
figurations considered in the dotted lines of Fig. [TOllead 
to the best fittings to the CMB multipoles, but the mat- 
ter powers of the same configurations (dotted line of the 
right panels) appear to be too large for any /c-value. For 
large spatial scales k < 0.05ft- Mpc^^, the observational 
data fit very well the predictions of the concordance 
model (solid lines); however, the predictions of the VT 
theory (dotted lines) are fully incompatible with the 
same data. The situation is much worse than in the 
concordance model. 

Best results have been obtained for erg > 0.9, r — 
0, and too large (small) values of Us {dns/d{ln fc)); 
namely, for unrealistic values of the above four param- 
eters. Let us now consider new values of Ofe/i^, which 
must be smaller than 0.0223. We have verified that (as 
it was expected) the Ce coefficients of three new config- 
urations similar to VTl, VTIO and VT13 but having 
Slft/i^ < 0.0223 are smaller than those of the proper 
VTl, VTIO and VT13 configurations. Then, start- 
ing from the new smaller coefficients, the same pro- 
cess -used above to look for the best fittings to the 
CMB and P{k) observed spectra- leads to as, r, Us 
and dns/d{ln k) values less realistic than in the case 
U,i)h? ~ 0.0223; hence, we cannot explain the observa- 



5 Evolution of vector perturbations in the VT 
theory 

In our version of the concordance model (GRl entry in 
Table [1]) we have assumed fife = —0.0012; nevertheless, 
this small curvature may be neglected. It follows from 
the effects produced by a much greater curvature with 
Ofc = —0.02. These effects (deviations of the dotted line 
with rerspect to the dotted-dashed one in the bottom 
panels of Fig. [9|) are small and, consequently, the ef- 
fects of a curvature with 17^ = —0.0012 are very small. 
On account of these comments and also for the sake of 
simplicity, we have taken U.^ — all along this section. 



Our flat model is chosen to be very similar to the 
concordance one. In order to ensure this similarity, all 
the parameters are assumed to be identical to those 
of the concordance model (see sectior(3l), excepting Ha 
and il^, whose values are assumed to be 0.7264 and 
0.001, respectively. This is the configuration VT16 of 
Table [1] In this case, it may be easily verified that 
the relation 57^ + ilc + + ~ 1 is satisfied, as 
it is mandatory in the flat case. Moreover, there is a 
dominant cosmological constant, and the dark energy 
due to is negligible. In spite of this fact, first order 
perturbations of A^ might produce important effects. 
In the absence of these perturbations the predictions of 
the configuration VT16 are very similar to those of the 
concordance model and, consequently, these predictions 
are compatible with observations. 

We are particularly interested in studying vec- 
tor perturbations, which could explain the WMAP 
anomalies in th e first CMB multipoles according to 
Morales &: Saej ( 20081) . without alterations of the 
power spectrum P{k) of the energy dens ity perturba- 
tions (scalar quantity, see BardeeiJ ( 1980l )). A study of 
these perturbations is performed in this section. In the 
linear regime, this study is independent of the existence 
of scalar or tensor perturbations. 

In VT theories, there are vector perturbations associ- 
ated to: the peculiar velocity Vi, the metric components 
hi — gcii, the vector compon ents A^, and the anisotropic 
stresses 11^ (|BardeeJl980h . As it is usually done, these 
stresses are assumed to be negligible. 

In a flat universe, vectors h, v, and A can be ex- 
panded in terms of the so-called fundamental vector 
harmonics, whose form is = exp(tfc ■ r ), wh ere k 
is the wavenumber vector (see iHu fc ^Yhitel (1199711) . A 
repre sentation of vectors e and e ' 
2007tl : 



( Morales k Saej|2007l) : (i) in the matter dominated 



(ifcifca/fc - ik2)/(jV2 , 



(10) 



era, 

Vc{T,k) = v^f^^^{k)/a{T) , 

B±(T,fc) = 6ilfo)f^m«,%)(fc)/fc'«'(r) (13) 

and, (ii) in the radiation dominated era, 

vf{T, k) = Vc{Teq, fc) — Constant , 

B±(T,fc) = 8prio)vtire,,k)/k^a^{T) , (14) 

where, r^q stands for the conformal time at matter- 
radiation equivalence. On account of Eqs.[T3landfT4l we 
conclude that quantities vf- are constant (decrease pro- 
portional to a~^) during the radiation (matter) domi- 
nated era. This is valid for any spatial scale. Let us now 
study the evolution of the gauge invariant quantities 
and A^ in the VT theory. 

The field equations of the theory couple the evo- 
lution of B^, v^, A^ , and the variables Aq and a 
describing the background. From the field equations 
of the VT theory, plus the perturbed line element 
ds'^ = a^i—dr'^ + Sijdx^dx^ — B^Qfdx''dT), and the 
perturbed vector field A^ = {Aq{t) , A^ Qf) , we have 
found the following equations describing -to first order- 
thc evolution of w^, B^, and A^: 



A^ 



fc2AoS± + (2- + 2^-fc2)A^ 



(15) 



2^B± = 
a 



'^[M2MB--A-) 
+Ao(Aoi3±-i±)] (16) 



(jMorales fc ^'^'^ 



^(A±-AoS±)] 



(17) 



= (ifcafcs/fc + ifci)/a\/2 (11) 



e± = ipa/fcv^ , (12) 

where a = {k\ + fc|)^/^. The expansions read as fol- 
lows h = B+Q+ + B-Q-, V = V+Q+ +v~Q-, and 
A = A+Q+ + A-Q-. Functions B^ik,T), v^{k,T), 
and A^(fc, r) describe the perturbation in momentum 
space. The differences v^{k,T) = ti^(fc, r) — B^{k,T) 
and the quantities A^(fc, r) are gauge invariant. 

In cosmological models based on GR (with and 
without cosmological constant), quantities A^ vanish 
and the time variations of and B^ are as follows 



From Eqs. ([3|) to ([8]) (after trivial modifications neces- 
sary to include vacuum energy) and (|15p to (|17p one can 
write a system of first order differential equations by us- 
ing appropriate variables and, then, this system can be 
numerically solved for suitable initial conditions. These 
conditions have been fixed in the radiation dominated 
era, at redshift z = 10*. Close to the initial redshift, 
we have proved that, approximately, all the variables 
involved in the problem evolve as powers of r. On ac- 
count of this fact, we have found the growing, constant 
and decaying modes, and we have chosen consistent ini- 
tial conditions. 

The evolution of quantity v'^ is given in Fig. Illl 
where dashed and solid lines correspond to GR and 
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VT, respectively. The unique difference between the 
two panels is the spatial scale L (k — 2tt/L). In the 
top panel, this scale is L = 2 x 10^ Mpc (superhorizon 
size) . We see that the separation between the two lines 
arises close to the end of the radiation dominated era. 
Moreover, the w^-values reaches the order 10~^^. After 
separation, the dashed (solid) line displays a decreas- 
ing (increasing) v'^ quantity. In the bottom panel, the 
spatial subhorizon scale is L = 2 x 10^ Mpc. In this 
case, it is evident that the solid line (VT theory) os- 
cillates around the dashed line (GR). We also see that 
the amplitude of the greatest oscillations is a few times 
10~^^, whereas the order 10~^^ is not reached. Oscil- 
lations (bottom panel) start later than the v'^ growing 
associated to the superhorizon scale (top panel). 

Fig. [13 shows the evolution of A'^. We see that 
this gauge invariant quantity decreases (oscillates) for 
superhorizon (subhorizon) scales. 

6 Discussion and conclusions 

In section [H the VT theory of section [2] has been ana- 
lyzed in detail for negligible ^'^-perturbations. Let us 
first summarize this analysis. 

There are pairs {^rn,h) explaining the SNe la ob- 
servations; however, CMB and matter power spectrum 
observations are not simultaneously explained for any 
of these pairs. 

Although we have studied a set of pairs covering the 
region of the {ilm,h) plane represented in Fig. IH in 
which SNe la observations are explained with p > 0.05, 
only the studies corresponding to three of these pairs 
have been described in detail. In these cases, main 
results are exhibited in Figs. ISltofTUl 

The main conclusions relative to the CMB angular 
power spectrum are now listed: (i) there is no a good fit- 
ting comparable to that of the concordance model in the 
full ^-interval (2,1150), (ii) the best fittings are shown 
in the dotted lines of Fig. [TOl These curves correspond 
to configurations VT9, VT12, and VT15, in which one 
has CTg > 0.9, T — 0, and too large (small) values of 
Us {dns/d{ln fc)); hence, these fittings are found for 
unrealistic values of the cosmological parameters, (iii) 
even in the case = 0.25 and h = 0.82 (the best 
fitting), the first observed multipoles {£ < 10) are not 
well explained, (iv) for other {flm,h) pairs, the fitting 
is worse, the number of unexplained small £ multipoles 
increases and the third peak is not well fitted; to see 
that, the middle left and bottom left panels of Fig. [10] 
must be compared with the top left one, which cor- 
responds to the best fitting, and finally (v) for the pa- 
rameter configuration VT9 leading to the best fitting to 



the CMB angular power spectrum, the SNc la data are 
only marginally explained with p ~ 0.07, and the Hub- 
ble constant is too large {h = 0.82). These conclusions 
strongly suggest that the VT theory under considera- 
tion does not work in the absence of perturbations. 

The study of the matter power spectrum also sug- 
gest that, in the absence of perturbations, the VT 
theory must be rejected. In fact, for the best fit- 
tings of the CMB angular power spectrum, the mat- 
ter power spectrum is not admissible for linear scales 
with k < 0.05ft- Mpc^^ (see the right panels of Fig. 
llOP : however, observations on these scales are very well 
explained by the concordance model without any bias 
between baryonic and dark matter (see Fig. [3]). 

According to the above comments, the VT theory 
appears to be inadmissible for any pair (r2m,/i) lying in 
the region represented in Fig. [6l provided that the en- 
ergy of the vector field A^^ plays the role of the dark en- 
ergy, there is no either vacuum energy or quintessence, 
and v4''-perturbations are negligible. Typical values of 
~ 1 — rim may be seen in Table [TJ 

A field whose contribution to the background en- 
ergy density of the universe is negligible does not sig- 
nificantly affect the background expansion and, con- 
sequently, it does not contribute to the SNe la dim- 
ming. However, the fluctuations of this field may 
produce crucial effects, for example, anisotropics in 
the CMB, which are fully absent in any homogeneous 
and isotropic background. These comments strongly 
suggest that fields being negligible at zero order -in 
cosmology- could be fully relevant at higher orders of 
perturbation theory. At first order, there are scalar, 
vector, and tensor modes and, if one or more of these 
perturbation components are not negligible, the field 
might be detected by the observation of the effects pro- 
duced by the non-negligible modes on the CMB and 
the matter power spectrum. This detection seems to 
be possible for rather small zero order energy densities 
Pj^. These arguments and those presented in section [T] 
-about the possible importance of the VT theories in 
the explanation of the CMB anomalies- have motivated 
the study presented in section [5| on the parameter con- 
figuration VT16 of Table [1] in which we have assumed a 
small background energy density with density parame- 
ter — 0.001. Thus, if both scalar and vector pertur- 
bations associated to A^^ are negligible, the small value 
of 17^ -in configuration VT16- leads to predictions in 
agreement with observations (see section [5]). 

After studying the VT theory in the absence of A^- 
perturbations, we have studied the evolution of vector 
perturbations in configuration VT16. This evolution is 
much simpler than that of the scalar perturbations and, 
moreover, both evolutions are independent (in the lin- 
ear regime); by this reason, we have been able to study 
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the evolution of the A'^ vector modes in section [5] The 
main conclusion is that, for superhorizon scales, quan- 
tities and grow. For the scale L = 2x 10'* Mpc, 
these quantities increase around 10^'^ orders of magni- 
tude from z — 10^ to present time; hence, even from 
small initial superhorizon perturbations, we can have 
large enough perturbations -during the recombination- 
decoupling process- producing anomalies in the small 
CMB £-multipoles. Subhorizon scales do not grow 
at the same rhythm and, furthermore, they oscillate; 
hence, only the superhorizon scales would be relevant 
at recombination-decoupling preventing anomalies for 
too large i values. The contributions of these super- 
horizon vector modes to the low i multipoles must be 
analyzed by using simulations and m ethods similar to 
those used bv [Morales fc Saed |2008l) . 

The scalar A'' -perturbations must be estimated in 
admissible inflationary models likely triggered by ap- 
propriate scalar fields coupled or not to field A^. On 
account of the results obtained in sections |4] and [5l 
any inflationary model leading to negligible scalar A^- 
perturbations must be discarded. The remaining mod- 
els might be admissible. It depends on the effects pro- 
duced by the A'^-perturbations after reheating. The 
evolution of the scalar perturbations is to be numer- 
ically studied with met hods similar to those used in 
standard cosmology (see iMa fc Bertschingeii (jl995l) ). 

Two are the main conclusions of this paper. The 
first one is important since it rules out, for the 
first time, some versions of the VT theory. In 
fact, we have concluded that the VT model studied 
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in previous papers (iBeltran Jimenez fc Marotd 12008 : 
Beltran Jimenez. Lazkoz fc Marotol bool " -where the 



A'^-energy density plays the role of the dark energy- 
is not compatible with current observations for negligi- 
ble A^'-perturbations. This conclusion must be taken 
into account in order to select inflationary models in 
the VT-theory. Our second conclusion is that model 
VT16 (see Table [1]) may explain the presence (ab- 
sence), at recombination-decoupling, of superhorizon 
(subhorizon) vector modes, which is the key to explain 
CMB anomalies in the low £ multipoles accor ding to 
the scheme proposed bv [Morales fc Saezj (|2008r ). It has 
been verified that there are growing vector modes in flat 
configurations different from VT16, which suggest that 
the existence of this kind of modes is a characteristic of 
the VT theories based on action ([ij . These conclusions 
motivate future investigations in the framework of the 
VT theory selected in section [21 and also in other VT 
theories. 
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Fig. 1 Solid line: supernova dista nce modulus, in terms of the redshift, z, for the concordance model. Observ ational data 
from the supernova legacy survey l|Astier et al.ll200^ ) and from a high redshift Hubble space telescope sample (|Riess et al.l 
l2006l ) are also shown. 




Fig. 2 Solid line: CMB angular power spectrum for the concordance model including lensing. Observational data (WMAP 
seven years) and error bars (including cosmic variance) are also shown. 
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Fig. 3 Solid line: all matter power spectrum estimated with CAMB for the concordance model. Observational data from 
the SDSS are also shown. 
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Fig. 4 Top: Power spectrum of the matter energy density perturbations in the concordance model. Solid line is the true 
power spectrum estimated with CAMB. Dashed line shows the power spectrum -for the same model- in the absence of 
baryon acoustic oscillations. Bottom: Ratio between the power spectra showed in the to panel 




Fig. 7 Quantity lOOp v.s. h for various values. Eacii curve corresponds to an Q^n value as it is displayed inside the 
panel. 
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Fig. 8 Left: angular power spectrum of the CMB temperature for some parameter configurations of Table[T] Right: power 
spectrum of matter fluctuations for the same cases as in the corresponding left panels. The configurations considered in 
each level are: GRl (solid), VTl (dotted), and VT2 (dashed), in the panels of the top level; GRl (solid), VTIO (dotted), 
and VTll(dashed) for the middle level; and GRl (solid), VT13 (dotted), and VT14 (dashed) for the bottom panels. 
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Fig. 9 Same as in Fig. [8] for the following configurations: GRl (solid), VT3 (dashed) and VT4 (dotted) in the top level, 
GRl (solid), VT4 (dotted), VT5 (dashed), and VT6 (dotted-dashed) in the panels of the middle level, and GRl (solid), 
VT6 (dotted-dashed) , VT7 (dotted), and VT8 (dashed) in the bottom panels. In the bottom right panel cases VT6, VT7 
and VT8 are indistinguishable. 
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Fig. 10 Same as in Fig. [8] for the following configurations: GRl (solid) and VT9 (dotted) in the top panels, GRl (solid) 
and VT12 (dotted) in the middle level, and GRl (solid) and VT15 (dotted) in the panels of the bottom level. In this figure, 
the same observational data as in Figs. [2] and [3] have been also displayed in the left and right panels, respectively. 
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Fig. 12 Same as in Fig. [TT]for tiie gauge invariant quantity A"*" 
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